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Abstract 

The author of this work have got explicit expressions for the timelike region at 
Next-to-Leading-Order (NLO) as for the coupling function so for QCD observables. 
These expressions were compared with approximate ones obtained in terms of "dou- 
ble logarithms". Next-Next-to-Leading-Order (NNLO) Pade approximation was also 
O \ discussed. 

In this paper some particular role of NLO was emphasised. It is related with the 
possibility to express any multi-loop expressions in terms of the two-loop ones. Thus 
t—i ', we can use NLO functions as a basis of some functional expansion, where higher loop 

' terms incorporation affects the expansion coefficients only. 

i ! 

1 Introduction 

cm : 

r** ■ Renorm-group equation (RGE) for QCD coupling function a s could be written as 

o ■ 

x-^- = -(P a 2 s +Pia 3 s +p 2 a 4 s + ...) (1) 
O dx 

with P n corresponding to (n+ l)-loop contribution. Its solution has a very simple form 

in the one-loop case: a^\x) = p \ nx - At NLO (two-loop approximation) it can be 

i | expressed through the Lambert function [|], ^, ^, [|] 

Ph. 

a^(x) = -^ - z = -HL x -i ( 2 ) 

^; ° s {x) Pii + w(z) ' Pie x ■ [ ) 

•rH , 

defined by a transcendental equation W(z)e w ^ = z. Here and further x = A is 
the scale parameter defining unphysical singularities positions. 

It should be mentioned that all RGE solutions obtained perturbatively have un- 
physical singularities such as "ghost poles". In the two- and three-loop cases the cut 
lying in infrared region emerges. Hence, the coupling behaviour in this region can't be 
described this way. Analytic Approach suggested four years ago by D.V. Shirkov and 
I.L. Solovtsov H ||, [7j allowed to obtain analytic expressions for the running coupling 
with the same ultraviolet asymptotic that ordinary RGE solutions have but with a 
stable infrared behaviour. Afterwards, this procedure was slightly corrected and gen- 
eralised by D.V. Shirkov |9|. In the Sec. 2 some basic facts from this approach were 
stated and also some two-loop calculations results were presented. 

In the Sec. 3 the Analytic Approach was applied for the two- and three- loop cases. 
It was proved in this work that we can obtain exact expressions for the two-loop case 
and Pade approximated third order /3-function. These results were compared with the 
ones obtained using "double logarithms" approximation. 

As follows from the publication [jl5|] it's possible to express the coupling constant 
for any-loop case as a power series of the two-loop coupling. It's also possible to 
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apply Analyticization Procedure to this expansion. The opportunity to express some 
arbitrary coupling function or observable in terms of two-loop functions is discussed in 
the Sec. 4. 



2 Modified Analyticization Procedure 

Standard Analytic Approach allows us to get rid of unphysical singularities and make 
it possible to research the coupling (and observables also) behaviour in the low energy 
region. Moreover it reduces the scheme dependence considerably. However it is also 
necessary to incorporate quark thresholds and to construct the procedure allowing to 
transfer the coupling (and observables also) to the time-like region. These problems 
were solved by D.V. Shirkov ||, |9| in his recent works. 



2.1 Q 2 -channel 

As it had been proposed earlier by the Analyticization procedure, it is possible to make 
the coupling analytic in a whole Q 2 complex plane except the cut on the real negative 
axis. However, this approach does not take into account the thresholds problem. An- 
alyticization procedure is based on Kallen-Lehmann spectral representation, it leads 
to the modified coupling function that is devoid of unphysical singularities. Besides, 
it has the correct asymptotic behaviour consistent with the perturbative result. The 
spectral representation has the form 

oo 


the spectral density can be taken as p(a) = Qa^ {— c) that is natural for the low 
energy processes analysis when there are three active quarks (/ = 3). Here, the new 
coupling function is defined through an imaginary part of usual, RG invariant, effective 
coupling a s continued on the physical cut. 

However, it is also possible to involve the threshold matching procedure. It is 
well known that one can relate the couplings for the regions with a different flavours 
numbers. One should match a{~ n ~ 1 (M 2 ) and a{~ n (M%) by the proper scale parameter 
redefinition. As a result, spectral density is discontinuous (8| and can be presented as 

p(a) = pf= 3 (a) + M 2 f )(p f (°) ~ P /= V))- (4) 

/>3 

Of course, the coupling obtained from this density is continuous. 



2.2 s- channel 

Some interesting experiments data correspond to the timelike region (s-channel) that 
leads to the necessity of the relevant redefinition of coupling function. Some approx- 
imate estimates, such as 2l(s) = a s (s) or 2l(s) = |a s (— s)|, were used earlier for this 
purpose. These estimates are suitable for qualitative analysis only in the region higher 
than 5GeV, and they can not reflect the correct analytical properties. The difference 
between these estimates and exact expressions is within several per cents in this region 
and this precision can be accurate enough for some calculations. 
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However, it is better to use the "dipole representation" for the Adler function in 
terms of an observable R(s) in the timelike region 

°«? 2 ) = e 2 / o ^»(.) (5) 
and an inverse expression 

R(s) = f ['^ ^dx. (6) 
2?r J_ s _ le x 

This transformation leads to the coupling function definition for s-channel. It can be 
implemented as 

2tt /_„_,, x v ; 



It can also be expressed in terms of the spectral density [10] 

-1 poo J 

2l( S ) = - / -p(a). (8) 
K Js cr 

These formulae have been applied to one-loop case by K. Milton and O. Solovtsova [11]. 
Corresponding expression for A(s) is 

1 1 In A 

2l = ^-r7rarctan-^ ( 9 ) 

Higher orders was also obtained in this paper. 



2.3 Thresholds 

As follows from the work by D.V. Shirkov ||, the effective s-channel coupling 2l(s) 
differs from 21-^ by a s/ii/it constant 

2l(s) = SI* + ^ at M| < s < M| +1 , (10) 

the s/iz/t constants can be easily calculated 

c/- 1 = c f + st^M)) - sa^-^Mj). (11) 

These constants values correspond to 4th, 5th and 6th quarks contributions. The shift 
constants provide the continuity of the analytic coupling function. Indeed, in spite 
of discontinuous spectral density we should get the continuous expressions from the 
spectral representation. 

Here are some typical values calculated at NLO: 

A 3 , MeV A4, MeV A 5 ,MeV A 6 ,MeV c 3 c 4 c 5 

300 251 174 72 0.0107 0.0028 0.00003 

400 346 249 108 0.0149 0.0035 0.00003 

Note that these contributions gather about 5 per cent to the effective coupling in the 
low energy region and can considerably influence the scale parameter value. 
Quite analogously, second shift constant and third shift constant are defined 

21; (s) = %{ + c{ at Mj < s < M| +1 , where (12) 
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c{- 1 = 4 + 2tf (M|) - *{-\Mj), i=2,3. (13) 
Some typical values of C2 and C3 are presented below: 

k A ,Mev c\ 4 4 4 4 4 

300 -0.0021 -0.0024 -0.0008 0.0026 0.0005 0.00002 
500 0.0007 -0.0024 -0.0004 0.0047 0.00009 0.00003 

3 Two-loop and three-loop applications 

3.1 Exact NLO results 

At NLO we can even obtain the analytical expressions for the s-channel Analyticized 
coupling function. As it was mentioned the two-loop renormgroup equation solution 
can be expressed in terms of the Lambert function (0). This representation also allowed 
us to obtain 2lj(s) functions exactly 

m = -^ + -z R -\ y (14) 

pi vr pia s [-s) 

2t2(s) = 4- Sln (l + 7T«s(- S )). (15) 
Trpi Po 



Ms) = 4-9 (Mi + § *.(-«)) - 

Pl 7T/3i [ Pq Pq 



(16) 



Ms) = (-£) 2 ^ {Mi + £*(-.)) - £*(-.) + |k(-)} , ^ (i 7) 

These expressions are obtained by using the contour integral (^) in the generalized 
form 

% n = / ^J-dx. (18) 



s—ie 



The integration can be implemented with rather simple variables substitutes. It's 
interesting that 2l n ( s ) can be presented as n — 2 residual terms of the Taylor expansion 
of — j^-Qdn(l + ^a s (— s)) over the powers of a s (— s) multiplied by — f^. 

3.2 Three- loop Pade approximation 

The renorm-group equation at NNLO can not be solved explicitly. However the Pade- 
approximated /3-function 

Ppade = -ml ( 1 + - A la 2 J = Poa 2 s + 0ia 3 s + /3 2 a* + ^a 5 s . . . (19) 
V ^-—aj Pi 

can be used in this case and the solution can be also obtained through Lambert func- 
tion g7| 



(3) / x Po 1 , 1 c -I h Pl P0P2 

l PadA x ) = -17--S — rr-— -, where z= -rex b , b = -3 , c = — ~ r 

- raae pil — c + W(z) eb /3q 0( 



Without Pade-approximation even in the three-loop case we can't get analytically 
the RGE solution for both coupling and observables. So it is natural to use the the 
expansion through two-loop functions to analyse multi-loop functions. 

One can obtain the s-channel coupling in this case exactly like in the two-loop case 

3.3 Using power expansion over the two-loop functions 
for NNLO analysis 

We can also get the three-loop functions without Pade approximation using the method 



suggested in the work [15|. The expansion ( J3T| ) allows us to obtain the following 
expressions 

% NNLO (s) = h^(s) + mf\s) + h&%\s) + fe 4 2ti 2) (s) + h*g\s) + ..., (22) 

for NNLO case k\ = 1, &2 = 0, k$ = k± = 0, k§ = This case differs from the 
case with Pade-transformed /^-function as the latter originates some "trail" of higher 
order terms. However this difference is quite small (when comparing new shift constants 
the difference between NLO and NNLO is about 15 per cents, while it is only 1 per 
cent between NNLO and "Pade" case). 



3.4 Comparing exact two-loop functions with approxi- 
mate ones 

Two-loop RGE solution (Q) is explicit and as it is argued in the present work the 
Lambert function is also very convenient for the timelike region research. In spite of 
these facts as a rule one takes some approximate expressions like 

<^> = sew { ( l ~ wW 1 + *fe -I) ' " i) } ' (23> 

As it was proved we can express this function in terms of @ 

o$ pp {x) = af\x) + k{x) (a^(x)) 4 + O [al{x)) (24) 
Here k(x) is some limited function. Some typical values of this function areQ: 

Q,MeV k(jjp) 

600 -0.7 
900 -0.78 
1200 -0.83 
1500 -0.85 
2100 -0.91 

This terms proportional to the fourth power of coupling function contributes for about 
4 per cent difference in lGeV region. We have to note that this can not be corrected 
by the scale parameter choice as it would originate nonzero coefficients before 2nd and 
3rd powers of in the expansion (S3) . 



1 We've taken A 3 = 300MeV 
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4 Obtaining the multi-loop functions in terms 
of two-loop ones 

4.1 Expansion 

Suppose some observable can be expressed as a power series of the coupling constant 

F{Q 2 ) = F.asiQ 2 ) + F 2 a 2 (Q 2 ) + F 3 a 3 s (Q 2 )... (25) 

However we face unphysical singularities, wrong analytic behaviour in the infrared (IR) 
region. We can not consider early attempts to express this observable in s-channel like 

F s (s) = F{a s (s) + F^a 2 s (s) + F& 3 s (s)... (26) 

to be very successful. Indeed, it involves the so called 7r 2 -terms and leads to the 
expansion coefficients augmenting [12, |l3|, this way also preserves unphysical sin- 



gularities. As a result of applying the Analyticization procedure the analytic coupling 
function and the functional expansion of observables both in Q 2 - and s- channels can 
be obtained 

F(Q 2 ) = FxAiQ 2 ) + F 2 A 2 (Q 2 ) + F 3 A 3 (Q 2 )... (27) 

F s (s) = Fi2l(s) + F 2 2l 2 (s) + F 3 2l 3 (s)... (28) 

After all, the Analytic approach leads to the scheme dependence reduction that is very 
important physical consequence. 

4.2 The expansion over two-loop functions 

As follows from the paper [15| we can obtain the expressions for multi-loop functions 
in terms of two- loop functions. Without analyticization we have the an expansion over 
the coupling function powers like 



..mult 



i-loop = a (2) + ^(2)2 + ^(2)3 + (29) 

Here we can easily apply the analytic approach. It leads to some modifications when 
we've got the functional expansion instead of the power one 

Amuiti-ioop = + k 2 A ( 2 2) + hA^ + ... (30) 

Thus two-loop functions can be considered as a minimal basis of any orders pertur- 
bation expansions. So any observable F (devoid of anomalous dimensions) can be 
presented over these functions in arbitrary loop order 

F = FxA® + F 2 A ( 2 2) + F 3 4 2) + .» (31) 

Some similar result (observables expansion as a power expansion over two-loop func- 
tions) was obtained independently by C. Maxwell [16|. 
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5 Conclusions 



In this paper the author has got an exact expressions for analytic running coupling 
at the two- loop order. These expressions differ from approximate ones by a term 
proportional to the fourth power of the coupling function. Timelike functions were 
also obtained just in the terms of Lambert function. As it was shown higher orders 
loop approximations can be obtained as an expansion over two-loop functions. We 
need to be only provided with /^-coefficients. 

This approach allows us to get obtain accurate theoretical predictions of the cou- 
pling function consistent with its analyticity properties. As it was argued the difference 
between exact expressions and approximate ones in the infrared region is high enough 
to be considered. 
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